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ABSTRACT 

The paucity of known High-Mass X-Ray Binaries (HMXB) with naked He donor stars (hereafter He 
star) in the Galaxy has been noted over the years as a surprising fact, given the significant number of 
Galactic HMXBs containing H-rich donors, which are expected to be their progenitors. This contrast 
has further sharpened in light of recent observations uncovering a preponderance of HMXBs hosting 
loosely bound Be donors orbiting neutron stars (NS), which would be expected to naturally evolve into 
He- HMXBs through dynamical mass transfer onto the NS and a common-envelope (CE) phase. Hence, 
reconciling the large population of Be-HMXBs with the observation of only one He-HMXB can help 
constrain the dynamics of CE physics. Here, we use detailed stellar structure and evolution models and 
show that binary mergers of HMXBs during CE events must be common in order to resolve the tension 
between these observed populations. We find that, quantitatively, this scenario remains consistent 
with the typically adopted energy parameterization of CE evolution, yielding expected populations 
which are not at odds with current observations. However, future observations which better constrain 
the underlying population of loosely bound O /B-NS binaries are likely to place significant constraints 
on the efficiency of CE ejection. 



1. INTRODUCTION 

As a star in a binary system evolves off of the main 
sequence (MS) and radially expands, it can overflow its 
Roche lobe and begin transferring mass onto its com- 
panion. If this transfer proceeds on a timescale shorter 
than the timescale in which the accretor can achieve ther- 
mal equilibrium, a common envelope (CE) develops and 
the stars begin to orbit through a combined atmosphere. 
Energy is then transferred from the binary orbit to the 
CE through fri ctional forces and torques which unbind 



the C E gases (Iben & Livio 1993 Taam & Sandquist 
2000 1 . This process has long Been discussed as the dom- 
mant mechanism for forming tight binaries from widely 
separated systems hosting one or two massive stars; it 
is thought to produce many cl osely interacting binaries 
such as Cata clysmic Variables (j Paczynski 19 761) and X 
Ra y Bina ries ( Eggleton fe Verbunt|19 86; Bailyn~& Grind 
lay|1987|. 



U nfortunately the hydrodynamics and long-term evo- 
lution during the CE phase are not well understood, mak- 
ing it difficult to accurately determine the efficiency of 
energy transfer from the binary orbit into the ejection of 
the CE. While advancements in computational hydrody- 
namics have elucidated several trends near the onset of 



the CE phase ( Rasio fe Livio 1996; Sandquist et al.|2000| 



Taam fc Ricker||2006[ |Ricker gjfaa m 2008), a simple 
parametrization must be employed in situations where a 
population of b inaries is to be considered. To this end, 
Webbink ( 1984 1 incorporated the myriad uncertainties of 
CE evolution into a single parameter acE, which governs 
the efficiency of transferring gravitational energy into the 
complete removal of the CE. Within this framework, a 
relation between the initial and final orbital separations 
can be written as: 



acE 



GM c M a G(M e + M c )M a 



= \E, 



bi n d | 



(i) 



where G is the gravitational constant, a$ and a/ are 
the initial and final orbital separations, M c , M e and M 
are the masses of the donor core, donor envelope, and 
accretor respectively, and ^bind is the energy necessary to 
unbind the CE. The value of ^bmd includes not only the 
(negative) gravitational binding energy, but also terms 
relating to the (positive) thermal energy of the plasma 
gas, t he ionization of H and He, and the disassociation 
of H 2 ( |Han et al.(l994| [l995| . 

Another uncertainty entering the calculation of ^und 



is the definition of the core-envelope boundaries ( Dewi 
fc Taurisj|2000| |Loveridge et al.||2010| ). As a rougrTcTeF 
inition of the stellar core, we assume the boundary to 
occur at the radius where the mass fraction of H drops 
below a critical value A min . Here we set X m ; n =0.1 and 
investigate the effect of changing this parameter in what 
follows. 

One long-standing expectation from models of binary 
evolution concerns the resilient population of high-mass 
X-ray binaries (HMXB) consisting of neutron stars (NS) 
accreting matter from the wind of naked He donors (He- 
HMXB). Such binaries are expected to form when the 
observed population of O/B - HMXBs evolves through a 
CE phase of the supergiant stellar component. The lack 
of observed He-HMXBs is particularly puzzling in light 
of recent observations uncovering an unexpectedly large 
population of HMXBs with Be-type stars (B-type stars 
which show emission- line (Be) spectra). Since each Be- 
HMXB contains a NS accretor along with a massive (08- 
B2) donor in a wide orbital period (P or b > 30 days), 
these systems are the expected progenitors of a bright 
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TABLE 1 



Name 


Spectroscopy 


cion \ \y / 


Pnvh (d) 


e 


7-Cas 


B0.5 IVe 


13.8 


203.59 


0.26 


0115+634 


BO. 2 Ve 


14.9 


24.3 


0.34 


0236+610 


B0.5 Ve 


13.8 


26.5 


0.55 


0331+530 


08-9 Ve 


23.0 


34.3 


0.3 


0352+309 


09.5 IIIe-BO Ve 


16.7 


250 


0.11 


0535+262 


BO Hie 


15.6 


111 


0.47 


0834-430 


BO-2 Ill-Ve 


12.0 


105.8 


0.12 


J1008-57 


09e-Ble 


15.6 


247.5 


0.66 


1417-624 


Bl Ve 


12.0 


42.1 


0.446 


J1946+274 


B0-1 IV- Ve 


13.8 


169.2 


0.33 


J1948+32 


B0 Ve 


15.6 


40.4 


0.03 


2030+375 


BOe 


15.6 


46.03 


0.41 


J2103. 5+4545 


B0 Ve 


15.6 


12.67 


0.40 



Table 1: List and parameters of the 13 galactic Bc-HMXBs with 
complete binary orbi tal information. Data obtained from Raguzova 
|fc Popov| j2005) and |Belczynski &c Ziolkowski| {2009;), and numerous 
references therein. The donor mass, orbital period and eccentricity are 
denoted by Md on , Porb, and e, respectively. Each donor mass is derived 
from the spectral classification using Tabic VIII of Habots & Hointzc 
41981} . 



HMXB population with He donors. Specifically, as both 
the orbital separation of observed Be-HMXBs is smaller 
than the star's supergiant radius and the mass ratio be- 
tween the Be donor and NS is large, Be-HMXBs in- 
evitably evolve into a CE. If the binary survives the 
CE event, the resultant system would contain a He star 
coupled with a NS in a relatively tight orbit. Since He 
stars ca n experience significa nt mass loss due to stellar 
winds (de Jager et al. 1988[ ), we would naively expect 
to observe such binaries as bright Hc-HMXBs. At least 
81 galactic Be-HM XBs are currently observe d , with 69 



reported by both Raguzova fe Popov (2005), Belczyn 
ski & Ziolkowskil (120091) , and numerous reference therein, 



as well as 12 additional s ystem s recently discovered by 



INTEGRAL flKaiir et al.||20 08 iRodriguez et al 



2010[|Corbet fc Krimm|2009p . Uf these systems, 13 have 



2009 



known orbital period, eccentricity, and spectral data, and 
are summarized in Table [T] Despite the large number of 
known Be-HMXBs, severaldecades of X-ray o bservations 



have uncovered only one He-HMXB, Cyg X-3 (|van Kerk- 
wijk et al.||1992 |. 



It is important to note that, throughout this paper, 
we specifically avoid the term Wolf-Rayet (WR) for the- 
oretically modeled He stars, as the exact definition of a 
WR star varies throughout the literature. While obser- 
vational studies focus on spectroscopic WRs (e.g. | Nugis| 
fe Lamers] 2000), many population synthesis investiga- 
tions treat all He st ars as WRs, applyi ng a WR-like 



wind mass-loss (e.g. Hurley et al. 2000 and Belczyn 
ski et aL"l|2008al who use the wind prescription from 



Hamann et al.||1995 Hamann fc Koesterke||l998| . On 
the other hand, other theoretical studies have set a min- 
imum mass for a WR in a HMXB of 5-8 M , consistent 



with the lowest WR mass inferred spectroscop i cally ( De 
Ponder et al. |1997 Vanbeveren et ah] 1998a Lommen 
et al.|2005p . The mismatch between the various defini- 



tions is further complicated by the fact that a low-mass 
He star (< 1.5M ) undergoing Roche-Lobe overflow onto 
a compact object may replicate the spectroscopic char- 
act eristics of the high er mass observed WR population 
( Lommen et al. 2005 1 . For this reason, in what follows 
we only call a star WR if it is a He star which shows spec- 
troscopic characteristics typical of observed WR stars. 



An analysis by van den Heuvel ( 1976 ) estimates a pop- 
ulation of He-HMXBs that exceeds the population of 
HMXBs with H-rich supergiant donors by a factor of 15, 
and argues that the lack of observed He-HMXBs provides 
an observational argument in favor of mergers during 
the dynamically unstable mass transfer phase. More re- 
cent p opulation synthesis calculations by |Lommen et aL] 
(2005) find that the galaxy is expected to contain ~1 
HMXB powered by the stellar wind of a He donor star 
larger than 7M , as well as ~1 HMXB powered by the 
Roche lobe overflow of a He star less massive than 1.5M Q . 
They conclude that it is equally probable that C yg X-3 is 
either type of system. Vanbeveren et al. ( 1998bh also pre- 
dict a significant population ot binaries hosting a NS and 
a He donor star, but note that the propeller effect may 
prevent significant accreti o n ont o the compact object. 



Van Bever & Vanbeveren (20001 note that the lack of 
HMXBs with He donors may be consistent with the lack 
of observed He sta rs in the solar neighborhood. Finally, 



Dewi et al. (12002) examined the production of double- 
1NS binaries from a population of low-mass He-HMXBs, 
finding that double-NS systems are only produced if the 
He star is undergoing He-shell burning during the pe- 
riod of RLO, while systems with He-core burning during 
RLO instead become white dwarf-NS binaries. However, 
they do not discuss the X-Ray characteristics of their 
He-HMXB population. 

Additionally, a significant effort has focused also on 
the population of observed runaway WRs, which may be 
expected to host compact object partners. |Moffat et al.| 
( |1982a|b | present spectroscopic observations of such sys- 
tems, but did not detect any X-ray bright source. Popu- 
lation synthesis studies targeting He stars more massive 
than 7 M found that only a small portion (1-2%) of the 
observed runaway binaries hosting an O-star and a He 
star may eventually evolve into He- HMXBs. Subsequent 
work by Vanbeveren et al. (1998a I found that at most 
3% of the galactic population ot massive He stars should 
be found in binaries containing a compact object. Both 
studies determined the paucity of observed He-HMXBs 
to be primarily due to the effect of disruptions due to 
the natal kick imparted to the compact objects at for- 
mation, as well as the probability of mergers during CE 
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events. However, these results cannot be directly applied 
to our study of He-HXMBs stemming from the observed 
Be-HMXB population, as we know that the system must 
survive as an intact binary through the NS natal kick. In 
any case, such population of runaway He stars is unlikely 
the evolutionary outcome of a population like the ob- 
served Be-HMXBs as the high spatial velocities of these 
systems are likely associated with the formation of the 
compact object and the observed O /B-HMXBs show sig- 
nificantly smaller spatial velocitie s ( Chevalier & Ilovaisky 
1998||van den Heuvel et al.||2000| ). 

In the present analysis, we use detailed stellar struc- 
ture and evolution models to investigate whether the 
discrepancy between Be-HMXB and He-HMXB obser- 
vations can be used to place constraints on the dynamics 
of the CE event, and specifically on oice- In § [2] we 
describe the modeling codes developed to calculate the 
orbital parameters of XRBs immediately before and af- 
ter the CE. In § [3j we use the currently observed sample 
of Be-HMXBs with measured orbital periods and eccen- 
tricities, and find that detailed calculations of massive- 
star binding energies and the typically adopted energy 
parametrization of CE evolution are consistent with the 
observed Galactic sample of one He-HMXB. In § |4j we 
investigate whether further limits can be placed on acE 
by simulating a grid of widely separated O /B-NS binaries 
and determining the survival probability of He-HMXB 
systems as a function of acE- In this way, we address 
the question of how a significant population of undis- 
covered O/B-NS binaries with wide orbits could further 
constrain energy deposition during the CE phase. We 
conclude in Section [5] with a discussion of how future ob- 
servations could be used to place further constraints on 
acE- 

2. SIMULATION AND MODELING CODES 

We have developed a detailed orbital evolution code 
suitable for investigating XRBs. This code tracks the 
evolution in time of the orbital separation and eccentric- 
ity of the binary, and the spin of the stellar component, 
accounting for the competing effects of stellar wind mass 
loss, wind accretion, tides exerted from the compact ob- 
ject onto the star, and angular momentum loss via grav- 
itational radiation. To account for changes in the stellar 
properties due to natural stellar evolution, the orbital 
evolution code is coupled to a detailed stellar structure 
and evolution code. The stellar evolution models are 
calculated with an up-to-date version of the stellar evo- 



lution code STARS | E ggleton| [19711 [19731 |Eggleton fc 
Kiseleva-Eggleton||20fe[|Pols et al = T|l995l ), while the INS 
is considered as a 1.4 M Q point mass. With this code, 
we follow the orbital evolution of Be-XRBs from the for- 
mation of the NS to the onset of Roche-lobe overflow 
(RLO), which corresponds to the beginning of the CE 
phase. At this stage we extract from the stellar model 
the donor's parameters (e.g. core and envelope masses, 
envelope binding energy, and spin) as well as the com- 
puted orbital separation and eccentricity. 
Tidal evolution is calculated in the standard weak- 



friction approximation ( Zahn||1977[ |1989[ ) , following the 
formalism of Hut (1981 ~ Specifically, we integrate nu- 
mericall y the set of differential equ ations as presented in 



For th is coefficient, we adopt a stellar model from |Claret| 
(20041 with a mass of ~ 16 M at solar metallicity and 
derive E 2 as: 



2.30 • tj^s 



logi; 2 = -5.71 - 2.51 • i™ 8 - 0.69 • t 7 ^ 

where £ms is time in units of the main sequence life- 
time. The fitting formula has only a very weak depen- 
dence on the initial mass in the observed range for the 
Be donor mass. The evolution of the orbital separation 
and ec centricity driven by s tellar wind are calculated fol- 
lowing Hurley et al. (2002). If the primary loses mass 



AMi, the orbit loses AMii?^ or b of angular momentum. 
If the secondary accretes some of this mass AM2, then 
AA^Rfilorb is returned to the orbit, where fl or b is the 
orbital frequency and Ri is the radius of the primary. 
The time evolution of the rotational freq uency driven by 



stellar wind follows Hurley et al.| ( 2000 ) and we assume 
that all the mass is lost uniformly at trie surface of the 
star. The evolution of the orbital separation and eccen- 
tric ity due to gravitational ra diation is calculated follow- 
ing Junker & Schaeferl ([ 1992 1). The accretion e fficiency is 
calculat ed according to Bondi & Hoyle ( 1944 ) following 
§ 4.2 of |Belczynski et al.| d2008b|). 



For each time step during trie orbital evolution we 
compute the Roche-lob e radius of the star at periastron 
(Sepinsky et al. 2007) and halt the calculation at the 
onset of the (JE phase. At this stage we extract the 
relevant binary and stellar parameters (see Eq. IT]) and 
calculate the resultant orbital parameters of the Dinary 
system from the energy prescription for CE evolution. If 
the radius of the donor star is found to exceed its Roche 
Lobe after the application of the CE prescriptions, we 
consider the system to end in a binary merger. 

Here we note that the mass-loss rate associated with 
the wind of a He star is uncertain and this, in turns, 
affects the X-Ray detectability of any surviving He-NS 
binary. Two very different models ar e quoted in the liter- 
ature. First, following the models of Hurley et al. (2000) 
we assign a single power law to the mass loss spectrum 
for all He stars stemming from massive binaries, given 
by: 



M, 



He 



1 x 1O- 13 (L/L ) 15 (3) 

However, work by Nugis fc Lamersj ( 2000 ) and |Nele-| 
mans fc van den Heuvel l |200l| ) haveproduced a broken 
power law where the mass loss declines prec ipitously for 
lower- mass He stars. The best fit given by |Dewi et al.| 
([2002] Eq. 6) follows: 



'2.8 x KT^L/L©) 1 - 5 log(L/L ) > 4.5 
~ S 4.0 x 10- 37 (L/£ Q ) 6 - 8 log(L/L Q ) < 4.5 



Mue. = 



3.3 of |Belczynski et "ah ( |2008b[ ), with the only mod- 



ification Denrgnii - flTfi - s^con^or^er tidal coefficient E2 



These models yield very different X-ray luminosities 
for the case of the lower mass Hc-HMXBs which are ex- 
pected to result from the observed Bc-HMXBs. In this 
work, we use both Eqs. ^ and Q and assess the effect 
on the observability of the He-HMXB populati on. After 
applying t hese w ind- mass loss rates, we follow [Belczyn-| 
ski et al.| ( |2008a[ Eqs. 39 & 83) to calculate the X-Ray 
luminosity of each individual system and to determine 
whether it enters the Chandra band. 
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Another uncertainty in this binary modeling concerns 
the evolutionary stage at which the Be-HMXB is ob- 
served, which affects how long the orbital evolution cal- 
culation should persist before the Be donor begins RLO. 
In order to constrain the error introduced by this uncer- 
tainty, we test three scenarios where current Be-HMXBs 
are assumed to be observed at ZAMS, terminal-age main 
sequence (TAMS), and an intermediate case in the mid- 
dle of the main sequence lifespan (MAMS). This effec- 
tively brackets the uncertainty stemming from the evolu- 
tionary state of the observed Be-HMXB population. We 
note that the long orbital periods of Be-HMXBs rule out 
the possibility that this MS evolution is interrupted by 
CE phases prior to the formation of the NS component, 
as CEs are known to produce binaries with substantially 
shorter orbital periods than observed in any Be-HMXB 
system ( Nebot Gomez-Moran et aT7||2011 ) . 

In addition, we investigate the production of He- 
HMXBs from currently observed Be-HMXBs by expand- 
ing on this procedure through the creation of a grid con- 
taining systems with parameters similar to the observed 
Be-HMXB population. Specifically, we use the detailed 
stellar evolution code STARS to create a grid of stellar 
models with masses between 1O-3OM (with a r esolution 
of lM^ ) and a probability distribution following | S alpet er 
(1955). We assume an initial orbital period in the range 
: 0-200d with a resolution of 1 day, and a density distribu- 
tion which is flat in the logarithm of orbital pe riod, which 
we use as a tracer for the orbital separation ( |Abt||1983| . 
In an alternative model we also investigate systems with 
orbital periods extending out to 400d and lOOOd, to de- 
termine the impact of this cutoff on our results. Finally, 
we employ a thermal eccentricity distribution following 
Heggie ( 1975 1. From this grid, the probability of a given 



progenitor surviving to become a visible He-HMXBs can 
be computed for an arbitrarily large population of likely 
progenitor systems. The normalization of the survivable 
probability to stellar environments can be ascertained 
through normalization against the observed number of 
Be-HMXBs as described in Section 21 

3. RESULTS FOR THE OBSERVED BE-HMXBS 

In Figure [T[ we illustrate the dynamics of our orbital 
evolution code by providing the detailed evolutionary 
history for a sing le simulation of 0236+610 (Table [If. 
Specifically, we plot the evolution in time of the or- 
bital separation, eccentricity, stellar radius, Roche-lobe 
radius, spin and orbital frequencies until the onset of 
RLO, when our calculation is halted. In this simulation 
we assume that the Be companion is currently observed 
at ZAMS. 

From this detailed analysis it is evident that the evo- 
lution of the each parameter is dominated by only the 
last 0.1% of its total lifetime, when the star develops 
a convective envelope which strongly enhance the tidal 
evolution. This suggests that simulations beginning at 
ZAMS, MAMS or TAMS will produce nearly identical 
binary parameters at the onset of the CE. 

The currently observed population of Be-HMXBs in- 
cludes only 13 systems with sufficient observational con- 
straints to allow for a detailed modeling of their orbital 
evolution forward in time (listed in Table [I]). Therefore, 
it is important to first determine which factors, other 
than CE dynamics, could prevent the observation of He- 
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FlG. 1. — Evolution of the orbital parameters for a binary with 
the properties of 0236+610 (see TableJTJ. Top: orbital separation 
(a), radius (R), and Roche-lobe radius (Rl)\ middle: spin and 
orbital angular frequency (u>); bottom: orbital eccentricity (e). The 
left panel shows the overall evolution, while the right panel shows 
a zoom-in at RLO. For the first ~ 13Myr, stellar wind mass-loss 
drives the evolution of the orbital separation, while the expansion of 
the star during its MS lifetime causes a decrease of the spin angular 
velocity. The evolution of the orbital eccentricity is driven by tides, 
which are too weak to affect it significantly. Towards the end of 
the MS the development of a convective envelope greatly increases 
the efficiency of tides. This mechanism controls the remaining 
orbital evolution leading to a decrease in the orbital separation and 
eccentricity, and driving the system into spin-orbit synchronism. 

HMXBs. Taking a small step aside, we test two likely 
factors, the relative lifetime of a bright He-HMXB phase 
compared to the Be-HMXB phase, and the expected X- 
Ray luminosity of the He-HMXB population. 

3.1. Alternative Explanations for the lack of Observed 
He-HMXBs 

One possible explanation for the lack of observed He- 
HMXBs concerns the relative duration of He-HMXB and 
Be-HMXB phases. If the time the stellar component 
spends as a He star is only a small fraction of the Be- 
HMXB lifetime, we would be unlikely to observe a large 
population of these systems regardless of the CE dynam- 
ics. In order to investigate this e ffect, we follow the cal- 
culations of Hurley et al. (20001 where the lifetimes of 
the MS and He-lVlS phases (tms an d Tffe respectively) 
are given by: 



tms = 



T~He 



1594 + 2707Af 4 + 146.6M 55 



0.04142Af 2 + 0.3426Af 7 
0.4129 + 18.81M^ e + 1.853M| 

We 



M 7 

-^Myr (5) 



■Myr (6) 



where Mo is the ZAMS mass and M# e is the mass of 
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— Full MS Age 

— MS Age Minus 5 Myr 

— MS Age Minus 10 Myr 
-Primary 3x as Massive 
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^ZAMS 

Fig. 2.— Relative lifetime of a 2.5 Mq (top left), 3 Mq (top right) 
and 4 Mq (bottom left) He star compared to the MS lifetime of O 
and B stars with various initial masses. To translate these results 
into a comparison of He-HMXB lifetimes and Be-HMXB lifetimes, 
we make four different assumptions regarding the fraction of the 
MS lifetime during which the star is a Be-HMXB: it's entire lifetime 
(black solid), beginning at 5 Myr (red dashed), beginning at 10 Myr 
(blue long dashed), at the SN age of a primary star which is three 
times as massive as M-zAMS (green dot-dashed). The mass range 
of 2.5-4 Mq for our population of He stars is set by the range of 
He core masses found in our models. 



the He main sequence star. We note that the direct com- 
parison of these lifetimes sets an extremely conservative 
lower bound on the population of expected He-HMXBs, 
as it assumes that the Be-HMXB is X-Ray bright for 
the entire MS lifetime of the donor star. This assump- 
tion is unrealistic, as it includes the portion of the MS 
lifetime which occurs prior to the evolution of the pri- 
mary star into a NS. In Figure [2] we plot the fractional 
lifetime for a He star of 2.5M Q , 3M©, and 4M© as a 
function of the ZAMS mass under four assumptions for 
the fraction of the total MS lifetime our observed sys- 
tems spend as a Be-HMXB. First, we assume that the 
system exists as a Be-HMXB for it's entire MS lifetime. 
Secondly, we subtract 5 Myr from the MS lifetime to 
account for the formation of an NS from the most mas- 



sive progenitors (see e.g. Muno et al. 2006). Third, we 
subtract 10 Myr from the MS lifetim e to account for th e 
average lifetime of NS progenitors (Hcger et al. 



Lastly, we follow the model of McSwam &z Gies 



2003 



2005 



who propose that the early evolution ot Be-HMXBs is 
governed by stable mass transfer from the NS progeni- 
tor onto the Be progenitor. This scenario sets an upper 
limit on the mass ratio between the primary and sec- 
ondary star at the onset of RLO based on the require- 
ment that the accreting star achieves thermal equilibrium 
on a timescale smaller than the mass transfer timescale 
of the primary star. While the exact mass ratio may 
depend sensitiv ely on the evolution ary state of each stel- 
lar compo nent ( Jben fc Livio|[l993 ) , for MS companions 
Hiellmin gJ |l989D set a range ot 2-4, while more recent 
work by Ivanova & Taam (2004) set a mass ratio of ap- 
proximately - ^ Motivated by these analyses, we lastly 
calculate the relative lifetime of the Be-HMXB phase as 
the MS lifetime of a Be star minus the MS lifetime of a 
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Fig. 3. — X-Ray Luminosity in the Chandra band for a HE-NS 
HMXB as a fu nction of the orbita l separation, using the He wind 
prescrip tion s of Hurley et al. (2000) (top) and the He wind prescrip- 
tions of |Dewi et al,| l |2UU^| (bottom) for He star masses spanning 
the range ot tnose produced via our stellar evolut ion code. In each 
case we use the X-Ray luminosity prescriptions of Bclczynski et al. 
< |2008a| 

progenitor star which is initially three times as massive. 
We note that all these scenarios are fairly conservative, 
due to the possibility that NSs in Be-HMXBs are formed 
via Electron-Capture supernovae (Nomoto 1984), which 
sets much stro nger constraints on the lifetime of the Be- 
HMXB phase (Linden et al. 20091. We note that a range 
of 2.5-4 M for the mass or the He core at the time of 
CE formation is strongly suggested by the results of our 
detailed stellar evolution models. 

Even considering the most massive He stars (with the 
shortest lifetime), as well as the most conservative calcu- 
lation of the Be-HMXB lifetime, we expect a fractional 
lifetime (t He/ T Be) of between 5-10% for Be stars be- 
tween 10-15 MQ. Thus the population of 81 currently 
observed Be-HMXBs predicts a population of at least 
6 He-HMXBs, which is at odds with the observation of 
only one such system at the 2a level. A significantly 
larger population exceeding 15 He-HMXBs is expected 
from more reasonable calculations of the luminous Be- 
HMXB lifetime and the He rich stellar mass. Thus, we 
may reject the hypothesis that the lack of observed He- 
HMXBs stems from their short lifespan. 

A second explanation for the lack of observed He- 
HMXBs concerns their assumed X-Ray luminosities. If 
mass lost from the He star is not efficiently transferred 
onto the NS, the systems may simply fall below the lu- 
minosity threshold of present observations. In this work, 
we assume that any system with an X-Ray luminosity 
in the Chandra band exceeding lxlO 34 erg s" 1 would 
have been detectable in galactic surveys (Andreas Zezas, 
Private Communication). We note that this assumption 
is conservative, and based primarily on the sensitivity 
of X-Ray survey missions such as ROSAT. A compari- 
son with the X-Ray luminosities of the Be-HMXB pop- 
ulation show sev eral systems with detecte d luminosities 
below this level (Raguzova fc Popov|[2005| . 

We note that these calculations depend sensitively on 



G 



both the mass of the He core determined by our evo- 
lutionary code, as well as the wind mass loss rate as- 
signed to these systems. We find that our He stars span 
a mass range of approximately 2.5-4 M with one out- 
lier exceeding 6 M Q . In Figure [3 we plot the calculated 
X-ray luminosity as a function of the orbital period for 
He-HMXBs hosting a He stellar component with a mass 
of 2.5M (black), 3M Q (red), and 4M Q (b lue), and fol- 
lowing the wind mass-loss m odels of both Hurley et al. 



( |2000| ) and |Dewi et"aLl ( |2002[ ) as giv en in Sec. jgj If 
winds following the prescription of IHurley ex al. 



stcll. 



ar 



winds following the prescription of Hurley et~~al. ( |2000 1 
are applied, these systems remain above the iummos- 
ity threshold out to orbital separations greatly exceeding 
those expected in post-CE binaries, which follow a log- 
arithmic normal distributio n spanning from 1.9h to 4.3 
days, with a peak at 10. 3h (Nebot Gomez-Moran et al. 
20111. However, in the case ot the weaker winds ap 



plied by Dewi et al. ( 2002 1 , the final orbital separation of 
the systems may affect its detectability in the Chandra 
band. Throughout what follows, we a dopt the more con - 



servative He mass- loss prescriptions of [Dewi et al. ( 2002 ) 
and only count as "detectable" those He-HMXBs with a 
luminosity exceeding lxlO 34 erg s _1 . 

Another process that might prevent He-NS binaries to 
be detected as X-ray so urces is the propeller effect . It has 



been previously noted (Vanbeveren et al. 1998b) that a 
high rotational velocity for the JN S accretor may prevent 
accretion of the donor wind material onto the NS. How- 
ever, this is unlikely to affect the He-HMXB population 
modeled here. Since these systems are observed to be X- 
Ray bright during the Be-HMXB phase, when the winds 
are less intense and the orbits are significantly wider, 
it is unlikely that the propeller effect will prevent the 
much stronger wind accretion during the tighter, post- 
CE phase. While the NS may be spun up during the CE 
phase, this spin-up is thought to be accompanied by a de- 
crease in the NS magnetic field. Such an effect is indeed 
observed for the pulsars contained in double-NS systems 
which have moved through CE evolution. These systems 
are expected to exist as the offspring of the He-HMXB 
population, and their reduced magnetic fields allow at 
least quasi-periodic accretion onto th e NS despite the 



high angular momentum in these NS (Romanova et al. 
2004 1 . In order to quantitatively examine the influence 
of the propeller effect on the He-HMXB population, a 
better understanding of NS spin and magnetic field evo- 
lution is required, but this is not possible at present. 

3.2. CE and the Observed Be-HMXB population 

Since the He-HMXB population has both a sufficiently 
long lifespan and high luminosity to indicate the exis- 
tence of numerous observable systems, another mecha- 
nism must halt the formation of these binaries. Since 
their progenitors are known to exist during the Be- 
HMXB phase, CEs stand as the only dynamical interac- 
tion which may eliminate He-HMXBs progenitors. Thus, 
it is possible to set constraints on the CE efficiency oice 
by demanding that enough of these systems are disrupted 
during a CE to bring the respective populations into 
agreement. Using the 13 Be-HMXBs with known orbital 
parameters as a template for the larger population of 
81 systems and employing the detailed stellar evolution 
models calculated with STARS, we employ our orbital 
evolution code to determine the orbital and stellar pa- 
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Fig. 4. — Number of surviving and observable (L^ > 10 34 erg s —1 ) 
He-HMXBs as a function of the CE parameter ciqe employing the 
orbital characteristics at RLO determined via our stellar evolution 
code for the 13 known Be-HMXBs. 

rameters of the observed Be-HMXBs at the onset of the 
CE. 

We again stress that we employ a calculation of |£"hmd| 
which includes not only the gravitational energy neces- 
sary to eject the primary envelope, but additionally the 
thermal energy in the envelope as well as H2 associa- 
tion and reionization. Since no other significant energy 
sources are available, the value of etcE must fall below 
unity in order to conserve energy. In Figure |4j we show 
the number of surviving He-HMXBs as a function of the 
CE efficiency ap e under the assumption that detectable 
Be-HMXBs were observed at the ZAMS, MAMS and 
TAMS, and with A min set to 15%, 10%, and 1%. We 
find nearly identical results in all scenarios, indicating 
that for o>ce ^ 1, mergers occur in all 13 of the observed 
Be-HMXB systems, regardless of their previous evolu- 
tionary history. Since one He-HMXB is in fact observed, 
and the lifetime of He-HMXBs may be smaller than that 
of Be-HMXBs, our model is only able to constrain acE 
to fall below unity, in line with theoretical expectations. 
We find that these results do indicate that current models 
are not missing any significant energy sources available 
to remove the CE, as the survivability of the CE phase 
exceeds 50% for olqe > 1-5, and these systems would be 
observable even w ith the weakened wind prescriptions of 
Dewi et all ((20021). 



4. A COMPLETE PARAMETER SPACE 

In the previous sections we have shown that CE-driven 
mergers are necessary to explain the discrepancy between 
Be-HMXB and He-HMXB observations and that the typ- 
ical CE parametrization, which demands an efficiency 
less than unity when all possible energy sources for un- 
binding the envelope are accounted for, is consistent with 
the observed sample of only one He-HMXB. However, 
these observed Be-HMXBs exist only as a subset for the 
potential class of He-HMXB progenitors. In addition to 
these luminous systems, there may exist a much larger 
population of binaries containing an O/B star and a NS 
in wide orbits, such that they are not bright X-ray sys- 
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Fig. 5. — Fraction of surviving and visibly detected He-HMXBs 
from a parameter space of initially widely separated O/B-NS bina- 
ries defined in the text with a maximum initial orbital separation 
between the CO and MS star of 200 days (solid black), 400 days 
(blue long-dash) and 1000 days (red short-dash). 

tems. This underlying population remains undetected 
because either the donor star does not carry enhanced 
winds stemming from the Be phenomenon, or the sys- 
tem is too widely separated for stellar material to be 
effectively accreted onto the NS. In any case, the binary 
dynamics of this underlying population are identical to 
the visible population, and evolution through a CE phase 
may similarly result in bright He-HMXBs. 

In order to model the evolution of these systems, we 
create a grid of binaries containing a NS and O/B donor 
following the parameters described in Section [2j In Fig- 
ure [5] we show the fraction of binaries which survives 
the CE yielding bright He-HMXBs. We assume max- 
imum initial orbital periods of 200d (black solid), 400d 
(blue dashed), and lOOOd (red dotted), following a distri- 
bution which is flat in the logarithm of the orbital period. 
We find that a potentially sizable fraction of He-HMXBs 
are created for larger values of acE, although the bounds 
depend strongly on the maximum assumed orbital period 
of the underlying O/B-NS population. We note the ob- 
servation of two Be-HMXBs with orbital periods above 
200d allows us to set this as an observed lower limit for 
our simulated population (see Table [7J) . We also note 
that we do not expect a significant variation in our re- 
sults if we were to vary the minimum orbital period in 
our sample grid (currently set at lOd). These relatively 
tightly bound systems are unlikely to survive a CE, and 
their inclusion in our models will not greatly affect the 
calculated number of He-HMXBs. 

In order to apply these results to the expected number 
of observed He-HMXBs, we must normalize the number 
of systems in our grid against the expected number of 
loosely bound O/B-NS binaries. We note that the num- 
ber of He-HMXBs expected from our simulated popula- 
tion of O/B-NS systems can be expressed as: 



Ntt«- 



He-HMXBs 



/ S (^)- X (— )No (7) 

N B T Be 

where i s is the survival probability of a given system 



from our simulation grid (shown in Figure pi), Nse/Ns 
provides the fraction of B-type donors which have Be- 
HMXB properties, Tff e / t Be—HMXB describes the rela- 
tive lifetime of the He-HMXB and Be-HMXB phases 
(shown in Figure [2]), and No is the observed number 
of Be-HMXB systems (No = 81 throughout this paper). 

The fraction of B-type stars which show emission-line 
(Be) spectra have been obse rved to vary between 2%- 
7% (McSwain & Gies 2005), although we note some 



sources have show n Be-fractions as high as 8.5% (Mo 
Swain et al.|2008| ). While this ratio may be substantially 



higher m binary systems if the spin-up of the Be popula- 
tion is due t o binary interactions, this line of reasoning is 
disputed by Oudmaijer & Parr (2010), who find a similar 
binary fraction tor both B and Be stars. In this work, 
we assume a Be-fraction of 7%, and a fractional lifetime 
for the He-HMXB population (THe/ T Be) of 20%, taking 
a central value from Figure [2] under the assumption that 
the primary progenitor was not more than three times 
as massive as the Be-star. We note that the X-Ray de- 
tectability of these systems is evaluated for each surviv- 
ing Hc-HMXB pro duced by our grid u sing the luminos- 

20021) and a luminosity 



ity prescriptions of |Dewi et al. 

cutoff of lxlO 34 erg s _i . From these values, we would 
anticipate a population of 230f, visible He-HMXBs. 

Thus, we constrain f s by comparing this expected pop- 
ulation of He-HMXBs to the observation of only a single 
system. Noting that a prediction exceeding four He- 
HMXBs would create a 2a discrepancy with observa- 
tion, we thus constrain the survivability of the CE phase 
to less than 2%. Comparison with Figure [5] thus con- 
strains the CE efficiency to a CE < 0.88, a CE < 0.75, 
and acE < 0.50 if the maximum orbital period is 200d, 
400d, and lOOOd, respectively. 

We note that the above calculation is conservative in 
several ways. First, we have assumed that all systems 
containing a Be star and NS are visibly bright X-Ray 
sources. Secondly, we have assumed that the observed 
population of Be-HMXBs can be translated to a popula- 
tion of O/B-NS with an orbital period which is logarith- 
mically flat starting at lOd. We note that the observed 
Be-HMXB population is instead biased towards systems 
with orbital periods around lOOd. While the lack of ob- 
served loosely bound systems may be due to luminosity 
cutoffs or simply to limited to observational time, the low 
period population is likely complete. This implies that 
the survival fraction of O/B-NS systems resembling the 
Be-HMXBs may be substantially higher. However, less 
conservative estimations are unlikely to significantly alter 
the constraint imposed on the CE evolution, as Figure [5] 
shows that the survivability of the CE phase plunges for 
smaller values of acE, implying that uncertainties in the 
estimation of the population of O /B-NS binaries has only 
a negligible effect on the number of expected He-HMXBs. 

5. DISCUSSION AND CONCLUSIONS 

Theoretical models predict the production of He- 
HMXBs through the CE evolution of widely separated 
binaries containing a NS and a massive donor. However, 
observations show a large population of Be-HMXBs and 
a significant lack of He-HMXBs. We find that detailed 
theoretical models predict the He-HMXB population to 
be sufficiently long-lived and luminous to be detected as 
the evolved offspring of the observed Be-HMXB popula- 
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tion. Noting that a CE phase acts as the only dynami- 
cal mechanism which may disrupt the production of He- 
HMXBs, we use these observations to set constraints on 
the CE efficiency parameter acE- Using the binary pa- 
rameters of the observed Be-HMXBs, we are only able 
to limit acE to fall below unity, echoing theoretical con- 
straints due to conservation of energy. Next, we simulate 
a larger grid of O /B-NS binaries with characteristics sim- 
ilar to the observed Be-HMXB population. From this 
grid, we constrain acE to be <0.88 for a population of 
O/B-NS binaries with a maximum orbital period of 200d, 
and possibly as low as uce <0.50 if the period extends 
to lOOOd. We note that this extension of the O/B-NS 
population to high orbital periods may have observa- 
tional support, as one X-Ray quiet system (B1259-63) 
with an orbital period of 1236d has been observed as a 
radi o pulsar with an optica l ly identified B-ty pe compan- 



ion (|Hughes fc Bailes|[1999| |Wex et al.|[1998 1 . 

We note that the stringency of our constraints are 
limited primarily by the low number of observed Be- 
HMXBs, and especially by the limited number of Be- 
HMXBs with known orbital period and eccentricity in- 
formation. We expect that observational detections of 
both new B-star NS binaries and determinations of the 
binary parameters of known Be-HMXBs will greatly re- 
duce these measurement errors and provide a more accu- 
rate understanding of CE evolution in massive binaries. 
Furthermore, we note that our findings are complemen- 
tary to previous studies which employ a combination of 
natal kicks and/or CE mergers to explain the small pop- 



Helium donors (|De Donder et al. 1997 


Vanbeveren et al. 


1998a 


Lommen et al.||2(J05|. Specihca 


liy, we extend the 



analysis to the lower mass range of He-HMXBs formed 



through the evolution of the Be-HMXB population and 
then use the number and binary properties of the ob- 
served Be-HMXBs to differentiate between the natal kick 
and CE merger hypotheses. We find that independent 
of any disruptions due to natal kicks (which would oc- 
cur prior to the Be-HMXB phase), CE mergers must be 
common in order to explain the relative paucity of He- 
HMXBs. Using this, we can place strict limits on the CE 
efficiency. 

Lastly, we note that a similar methodology may be 
applied to the population of known NS-NS binaries con- 
taining a pulsar, in order to determine whether the rate 
of NS-NS production is itself consistent with the low 
survivability probability assigned to Be-HMXBs moving 
through CE phases. Since these systems would addi- 
tionally experience a He-HMXB phase in between the 
Bc-HMXB and NS-NS phases, the existence of a large 
NS-NS population inconsistent with the small number of 
Be-HMXBs would instead point towards the existence of 
an X-Ray quiet population of He-HMXBs, and may be 
used as a further test of the results obtained here. 
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